Linear Relationship Statistics in Diffusion Limited Aggregation 
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We show that various surface parameters in two-dimensional diffusion limited aggregation (DLA) 
grow linearly with the number of particles. We find the ratio of the average length of the perimeter 
and the accessible perimeter of a DLA cluster together with its external perimeters to the cluster 
size, and define a microscopic schematic procedure for attachment of an incident new particle to the 
cluster. We measure the fractal dimension of the red sites (i.e., the sites upon cutting each of them 
splits the cluster) equal to that of the DLA cluster. It is also shown that the average number of the 
dead sites and the average number of the red sites have linear relationships with the cluster size. 

PACS numbers: 61.43.Hv, 64.60.al, 68.35.Fx, 47.57.cb, 61.05.-a 



Diffusion-limited aggregation (DLA) is a model of 
a growing cluster, originally proposed by Witten and 
Sanders [1] . The model has been shown to underlie many 
pattern forming processes including dielectric breakdown 
[2], electrochemical deposition [3], viscous fingering and 
Laplacian flow [4] etc. It is defined by a simple stochastic 
model on a square lattice as follows. A seed particle is 
located at the center of the lattice, and then a random 
walker is released from infinity — operationally, from a 
point of radius much larger than the radius of the grow- 
ing cluster. Upon contacting, the random walker sticks 
irreversibly to the cluster. Repeating the process, leads 
to an intricate ramified structure whose surface in the 
plane grows proportional to the bulk (this will be shown 
in this letter). 

This procedure is equivalent to solving Laplace's equation 
outside the aggregated cluster with appropriate bound- 
ary conditions. In two dimensions, since analytic func- 
tions automatically obey Laplace's equation, the theory 
of conformal mappings provides another mechanism for 
producing the shapes. This method has been directly 
used by Hastings and Levitov to study DLA [5] . 
One of the most interesting aspects of such an aggregate 
is the multifractal behavior of the growth site probability 
distribution (the harmonic measure) {pi}, where is the 
probability that the site i, belonging to the perimeter of 
the cluster, will grow at the next time [6, 7] . The screened 
sites with tiny growth probability, play an important role 
in determining the multifractality, while evaluating the 
harmonic measure on these sites is too difficult. The the- 
oretical difficulty emerges from solving equations with 
boundary conditions on a complicated growing interface. 

Different numerical methods applied for large aggre- 
gates in plane, show that DLA does not follow a simple 
fractal pattern and deviates from linear self-similarity [8] . 
However most of the studies have been focused on the 
scaling behavior of different quantities in DLA, here in 
this letter, we show that there exist couple of linear re- 
lationships which, to our knowledge, have not been ad- 



dressed yet. We show that the perimeter and the accessi- 
ble perimeter of a growing aggregate, together with their 
externals grow linearly with the cluster size n, with dif- 
ferent rates. The number of sites which are not accessible 
for the incident random walker will be shown to increase 
linearly with size. Our results indicate that the almost 
particle attachments (which a particle is assumed to be at 
the center of a square of the mesh size a), choose single- 
side contact. We also use the concept of red sites, as a 
measure of ramification of the aggregates and compute 
their fractal dimension. The average number of the red 
sites is another quantity which has a linear relationship 
with the number of aggregated particles. 

In order to investigate the behavior of the aforemen- 
tioned quantities, we simulated several independent on- 
lattice DLAs of different mass up to n = 10^ parti- 
cles. One can see that an ensemble of simulated DLAs 
is strongly fiuctuating, and to obtain consistent results 




FIG. 1: (Color online) A putative small DLA cluster of n 
aggregated particles shown in shaded squares. The marked 
squares with x show the dead sites which are inaccessible for 
the random walker coming from infinity. The hull of the DLA, 
defined in the text, is the union of the solid and dashed lines. 
The outer perimeter is the union of the solid and dotted lines. 
The accessible perimeter is the hull of the union of the DLA 
cluster and the dead sites. The (n -I- l)-th incident random 
walker may have either positive, negative or zero contribution 
to the length of the hull, denoted hj p+, p_ or po, respectively. 
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the averages have to be taken over a large number of 
samples. In this letter, the averages are taken over 5000 
independent samples for different cluster sizes. 

The first quantity that we measure is the perimeter 
or hull of an aggregate. Consider a cluster of aggre- 
gated particles on a square lattice (see Fig. 1), where 
each frozen particle in the cluster is assumed to be lo- 
cated at the center of a plaquettc. To define the hull, a 
walker moves clockwise around the aggregate and along 
the edges of the corresponding lattice starting from a 
given boundary edge of the cluster. The direction at 
each step is always chosen such that walking on the se- 
lected edge leaves a frozen particle at the right and an 
empty plaquette at the left of the walker. If there are two 
possible ways, the preferable direction is one that is on 
the right of the walker. The directions right and left are 
defined locally according to the orientation of the walker. 
This algorithm yields the hull of the DLA whose length I 
(in the units of the mesh size a), is equal to the number 
of steps until the walker returns to the starting edge. 
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FIG. 2: (Color online) The hull (dark colored) and the accessi- 
ble perimeter (light colored) of a DLA cluster of size n = 10^ . 



Let us denote the length of the hull of a cluster of size 
n, by We now show that In must be proportional to 
the cluster size. 

To have a cluster of size n-f-l, a random walker is re- 
leased from infinity. It sticks to an edge of the hull 
with the probability proportional to the harmonic mea- 
sure there. We define the three possibilities as shown 
in Fig. 1, during which the random walker can stick to 
the cluster having either single^ double or triple-side con- 
tacts. Upon selecting each of them, it may have positive, 
negative or zero contribution to the length of the up- 
dated hull. According to its contribution, we denote the 



probability that the random walker selects each of the 
three possibilities as p+, p_ or respectively. It can be 
easily checked that the events with single side contacts 
have always a positive contribution of -1-2 to the length 
of the hull, while the two other possible events can either 
have zero or negative contribution of 2n_, with ?7_ > 1. 
One can indeed define the length of the updated hull by 
considering these contributions to the previous length 
using the following recursive relation 



In+l = In + 2[p+ - 



(1) 



As our following experiments suggest, for n ^ 1, the last 
term in the above relation seems to be independent of 
the cluster size and therefore, one can obtain that In ~ 
2[p-|_ — n^pJ\n. Our simulation result for the average 
length of the hull Z„ as a function of the cluster size n is 
in a good agreement with this linear relationship, see Fig. 
3. We find that ^ = 1.791(2). This yields the infimum of 
the probability that the incident random walker chooses 
a single side contact, and hence we obtain > 90%. 

In addition to studying the hull of the aggregates, we 
also study the external perimeter of the hull which can 
measure the number of sites got trapped in the fjords 
(proportional to n_). To define the external perimeter, 
we first close off all the narrow passageways of a lattice 
spacing on the DLA cluster and then looking at the hull of 
the resulting cluster (see Fig. 1). We find experimentally 
that the length of the external perimeter l^ has also a 
linear relationship with the cluster size n. The best fit 

to our data shows that ^ = 1.325(2). 

All sites on the hitherto considered perimeters are not 
necessarily accessible by the incident random walker in 
all regions of the aggregate. Therefore, it is of interest 
to measure the totally accessible perimeter which is, in 
principle, a hull that surrounds the union of the DLA 
cluster and all inaccessible sites by the incident random 
walker coming from infinity— see Fig. 1. 
In order to determine the accessible perimeter, we intro- 
duce an algorithm which seems to be more efficient than 
one used in ref. [10]. This may be called '"burning al- 
gorithm\ during which each accessible site around the 
cluster will be marked as a burning site, and all cluster 
sites and inaccessible ones will be left unmarked. 

The algorithm begins by drawing a box that includes 
the entire cluster without touching it, and 'marking in' 
from the boundary of the box. A boundary site of the 
box is selected and is marked as '/ree' if there is not any 
cluster site at its nearest neighborhood and as ^burnV 
otherwise. If the site left 'free '-marked, after checking all 
of its nearest neighbor sites and marking each of them 
as 'free' or 'burnt' as before, it is recolored as 'burnt'. 
Repeating this procedure for all 'free'-marked sites and 
all unmarked boundary sites of the box, partitions inside 
the box into two regions, one containing the sites which 
are being marked as 'burnt' sites and are all accessible by 
the random walker, and the other region i.e., inaccessi- 
ble region, is the union of the cluster sites and unmarked 
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sites (or dead sites) which are not accessible by the ran- 
dom walker. The accessible perimeter is then the hull of 
the inaccessible region. 

We find that the average length of the accessible perime- 
ter and the average number of dead sites iV^, grow lin- 
early with the cluster size n. The linear relations can be 

obtained as shown in Fig. 3. according to ^ = 1.203(2) 
and = 0.554(2), respectively. 

Following the same ratiocination as Eq. (1) for these 
observations, one can estimate the probability that the 
incident particle sticks to the accessible perimeter of DLA 
with a single-contact. We denote the same probabilities 
as before that the incident random walker coming from 
infinity to stick to the accessible perimeter with single, 
double or triple- side contacts by p^, p'_ or p'q, respec- 
tively. Note that for the hull of accessible perimeter, all 
single, double or triple-side contacts always have posi- 
tive, zero or negative contribution (of +2, or —2, i.e., 
n_ = 1) to the length, respectively. Unlike DLA perime- 
ter (because of screening with the dead sites) , all sites on 
the accessible perimeter have nonzero growth probability. 
It can be again shown that ~ 2[p^ — p'_]n. This re- 
sult is checked in Fig. 3, according which we obtain that 
the length of the accessible perimeter ^" grows linearly 

with the cluster size n as ^ = 1.203(2). It can thus be 
inferred that the incident square-like random walker at- 
taches to the active zone of the cluster with probability of 
p'+ ^ 60% by choosing the single-side contact. Since all 
fjords of narrow throat are filled by the dead sites inside 
the accessible perimeter, so we can estimate the number 
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FIG. 3: (Color online) From top to bottom: the rescaled 
length of the perimeter external perimeter l„, accessible 
perimeter l„, accessible external perimeter , and rescaled 
number of the dead sites Nd versus the cluster size n. The 
errors are less than the symbol size. The solid lines indicate 
the best linear fit of zero slope. 
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FIG. 4: (Color online) Linear relationships between different 
perimeter lengths and the area enclosed by them. From top 
to bottom: the length of the external perimeter l„, accessible 
perimeter l„ and accessible external perimeter Z„ , rescaled by 
their enclosed area denoted by A„, and , respectively. 



of sites A/L on it with triple-side contact possibility, by 
closing off all the narrow passageways of a lattice spacing. 
This yields the outer of the accessible perimeter whose 
length Z"' has also linear relationship with the size of he 

cluster, i.e., as shown in Fig. 3, ^ = 1.068(2). Com- 
paring with the same relation for I", one can conclude 
that M- ^ 0.067n. However, by closing off the narrow 
passageways, all possible remaining contacts will be sin- 
gle or double-side ones, nevertheless it is not possible to 
estimate an exact value for It is because that these 
attachments can also have negative contribution to the 
length of the external accessible perimeter. 

These linear relationships seem to be the characteris- 
tic features of DLA clusters. In fact, for common fractals 
appear in two-dimensional statistical mechanics, such as 
critical Ising or percolation clusters, the length of the 
cluster boundaries (or loops) I with the fractal dimen- 
sion of c?/, has a scaling relation with the area of the 
loops as ^ P^^'^f . 

Motivated by this relation, we examine the behavior of 
the length of the perimeters versus the area enclosed by 
them. As depicted in Fig. 4, we find that all the perime- 
ter lengths Z„, Z„ and /„ , have linear relations with their 
area. Note that this area is not necessarily equal to the 
cluster size (multiplied by the square lattice spacing a^). 
These linear relationships show that applying the scal- 
ing relation A ~ p/'^f to measure the fractal dimension 
of the perimeter of DLA clusters would thus lead to a 
misleading result reported in [9]. 

The fractal dimension of the perimeters df can be com- 
puted by using the scaling relation between the average 
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FIG. 5: (Color online) Main frame: the average number of 
red sites A^r versus the gyration radius Rg of DLA clusters 
of different sizes. Inset; linear relation between Nr and the 
cluster size n. 



length of the perimeter I and a linear size scale, e.g., 
the gyration radius Rg, i.e., / ~ Rg ■ As long as the 



gyration radius of the growing cluster is used as the lin- 
ear size scale in this relation, we find that the fractal 
dimension of all perimeter lengths, within the statisti- 
cal errors, are equal with the same value as the frac- 
tal dimension of DLA cluster. This result is in agree- 
ment with the same one reported in [10]. Nevertheless, 
if we scale the perimeter lengths with their gyration ra- 
dius (i.e., the gyration radius of the loops produced by 
the perimeters themselves), we find a minor difference 
among the fractal dimensions. If dy-, d/, d^, d^ and dy 
represent the fractal dimension of the DLA cluster, the 



perimeter, the external perimeter, the accessible perime- 
ter and the accessible external perimeter, respectively, we 
find that d^ = 1.707(3), d/ = 1.710(3), d^ = 1.717(3), 

d'^ = 1.723(3) and dj 1.725(3). 

In the following, we use the concept of the red sites 
borrowed from two-dimensional critical structures [11], 
to have a quantitative measure for the ramification of 
DLA cluster. A red site on the DLA cluster denotes a site 
that upon cutting leads to a splitting of the cluster. This 
measures the number of nodes connected by effectively 
one-dimensional links. 

We carried out simulations in order to compute the 
fractal dimension of the red sites d,, on DLA cluster by 
using the scaling relation Nr ~ Rg'', where Nr is the 
number of red sites and Rg is the gyration radius of the 
DLA cluster. Due to the large amount of time, the simu- 
lation was performed for clusters of size n < 4.5 x 10"*. As 
shown in Fig. 5, we find that, within statistical errors, 
the fractal dimension of the red sites is equal to that of 
the DLA cluster, i.e., dr = 1.709(3). 
We also find that the number of the red sites grows 
linearly with the cluster size according to the relation 
^ = 0.673(2). 

In conclusion, we found that different surface parame- 
ters on DLA cluster such as the length of the hull, the ac- 
cessible perimeter and their externals and also the num- 
ber of dead sites grow linearly with the cluster size. These 
have been used to investigate the microscopic feature of 
the cluster growth by measuring the probability that a 
square-like incident random walker attaches to the clus- 
ter by choosing either the single, double or triple-side 
contacts. We also found that the border of DLA grows 
linearly with the total area enclosed by it. We have mea- 
sured the fractal dimension of the red sites on the DLA 
cluster and found it to be equal with that of the clus- 
ter itself. The average number of the red sites has been 
shown to have a linear relationship with the cluster size. 
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